Split Pool Ligation-based Single-cell Transcriptome sequencing (SPLiT-seq) data processing pipeline comparison.

BACKGROUND: Single-cell sequencing techniques are revolutionizing every field of biology by providing the ability to measure the abundance of biological molecules at a single-cell resolution. Although single-cell sequencing approaches have been developed for several molecular modalities, single-cell transcriptome sequencing is the most prevalent and widely applied technique. SPLiT-seq (split-pool ligation-based transcriptome sequencing) is one of these single-cell transcriptome techniques that applies a unique combinatorial-barcoding approach by splitting and pooling cells into multi-well plates containing barcodes. This unique approach required the development of dedicated computational tools to preprocess the data and extract the count matrices. Here we compare eight bioinformatic pipelines (alevin-fry splitp, LR-splitpipe, SCSit, splitpipe, splitpipeline, SPLiTseq-demultiplex, STARsolo and zUMI) that have been developed to process SPLiT-seq data. We provide an overview of the tools, their computational performance, functionality and impact on downstream processing of the single-cell data, which vary greatly depending on the tool used. RESULTS: We show that STARsolo, splitpipe and alevin-fry splitp can all handle large amount of data within reasonable time. In contrast, the other five pipelines are slow when handling large datasets. When using smaller dataset, cell barcode results are similar with the exception of SPLiTseq-demultiplex and splitpipeline. LR-splitpipe that is originally designed for processing long-read sequencing data is the slowest of all pipelines. Alevin-fry produced different down-stream results that are difficult to interpret. STARsolo functions nearly identical to splitpipe and produce results that are highly similar to each other. However, STARsolo lacks the function to collapse random hexamer reads for which some additional coding is required. CONCLUSION: Our comprehensive comparative analysis aids users in selecting the most suitable analysis tool for efficient SPLiT-seq data processing, while also detailing the specific prerequisites for each of these pipelines. From the available pipelines, we recommend splitpipe or STARSolo for SPLiT-seq data analysis.

Species-specific responses during Seoul orthohantavirus infection in human and rat lung microvascular endothelial cells.

Seoul orthohantavirus (SEOV) is a rat-borne zoonotic virus that is transmitted via inhalation of aerosolized infectious excreta, and can cause hemorrhagic fever with renal syndrome (HFRS) in humans worldwide. In rats, SEOV predominantly exists as a persistent infection in the absence of overt clinical signs. Lack of disease in rats is attributed to downregulation of pro-inflammatory and upregulation of regulatory host responses. As lung microvascular endothelial cells (LMECs) represent a primary target of infection in both human and rats, infections in these cells provide a unique opportunity to study the central role of LMECs in the dichotomy between pathogenicity in both species. In this study, host responses to SEOV infection in primary human and rat LMECs were directly compared on a transcriptional level. As infection of rat LMECs was more efficient than human LMECs, the majority of anti-viral defense responses were observed earlier in rat LMECs. Most prominently, SEOV-induced processes in both species included responses to cytokine stimulus, negative regulation of innate immune responses, responses to type I and II interferons, regulation of pattern recognition receptor signaling and MHC-I signaling. However, over time, in the rat LMECs, responses shifted from an anti-viral state towards a more immunotolerant state displayed by a PD-L1, B2M-, JAK2-focused interaction network aiding in negative regulation of cytotoxic CD8-positive T cell activation. This suggests a novel mechanism by which species-specific orthohantavirus-induced endothelium and T cell crosstalk may play a crucial role in the development of acute disease in humans and persistence in rodents.

Distinctive cell-free DNA methylation characterizes presymptomatic genetic frontotemporal dementia.

OBJECTIVE: Methylation of plasma cell-free DNA (cfDNA) has potential as a marker of brain damage in neurodegenerative diseases such as frontotemporal dementia (FTD). Here, we study methylation of cfDNA in presymptomatic and symptomatic carriers of genetic FTD pathogenic variants, next to healthy controls. METHODS: cfDNA was isolated from cross-sectional plasma of 10 presymptomatic carriers (4 C9orf72, 4 GRN, and 2 MAPT), 10 symptomatic carriers (4 C9orf72, 4 GRN, and 2 MAPT), and 9 healthy controls. Genome-wide methylation of cfDNA was determined using a high-resolution sequencing technique (MeD-seq). Cumulative scores based on the identified differentially methylated regions (DMRs) were estimated for presymptomatic carriers (vs. controls and symptomatic carriers), and reevaluated in a validation cohort (8 presymptomatic: 3 C9orf72, 3 GRN, and 2 MAPT; 26 symptomatic: 7 C9orf72, 6 GRN, 12 MAPT, and 1 TARDBP; 13 noncarriers from genetic FTD families). RESULTS: Presymptomatic carriers showed a distinctive methylation profile compared to healthy controls and symptomatic carriers. Cumulative DMR scores in presymptomatic carriers enabled to significantly differentiate presymptomatic carriers from healthy controls (p < 0.001) and symptomatic carriers (p < 0.001). In the validation cohort, these scores differentiated presymptomatic carriers from symptomatic carriers (p ≤ 0.007) only. Transcription-start-site methylation in presymptomatic carriers, generally associated with gene downregulation, was enriched for genes involved in ubiquitin-dependent processes, while gene body methylation, generally associated with gene upregulation, was enriched for genes involved in neuronal cell processes. INTERPRETATION: A distinctive methylation profile of cfDNA characterizes the presymptomatic stage of genetic FTD, and could reflect neuronal death in this stage.

Runx1+ vascular smooth muscle cells are essential for hematopoietic stem and progenitor cell development in vivo.

Hematopoietic stem cells (HSCs) produce all essential cellular components of the blood. Stromal cell lines supporting HSCs follow a vascular smooth muscle cell (vSMC) differentiation pathway, suggesting that some hematopoiesis-supporting cells originate from vSMC precursors. These pericyte-like precursors were recently identified in the aorta-gonad-mesonephros (AGM) region; however, their role in the hematopoietic development in vivo remains unknown. Here, we identify a subpopulation of NG2+Runx1+ perivascular cells that display a sclerotome-derived vSMC transcriptomic profile. We show that deleting Runx1 in NG2+ cells impairs the hematopoietic development in vivo and causes transcriptional changes in pericytes/vSMCs, endothelial cells and hematopoietic cells in the murine AGM. Importantly, this deletion leads also to a significant reduction of HSC reconstitution potential in the bone marrow in vivo. This defect is developmental, as NG2+Runx1+ cells were not detected in the adult bone marrow, demonstrating the existence of a specialised pericyte population in the HSC-generating niche, unique to the embryo.

Comparison of Single Cell Transcriptome Sequencing Methods: Of Mice and Men.

Single cell RNAseq has been a big leap in many areas of biology. Rather than investigating gene expression on a whole organism level, this technology enables scientists to get a detailed look at rare single cells or within their cell population of interest. The field is growing, and many new methods appear each year. We compared methods utilized in our core facility: Smart-seq3, PlexWell, FLASH-seq, VASA-seq, SORT-seq, 10X, Evercode, and HIVE. We characterized the equipment requirements for each method. We evaluated the performances of these methods based on detected features, transcriptome diversity, mitochondrial RNA abundance and multiplets, among others and benchmarked them against bulk RNA sequencing. Here, we show that bulk transcriptome detects more unique transcripts than any single cell method. While most methods are comparable in many regards, FLASH-seq and VASA-seq yielded the best metrics, e.g., in number of features. If no equipment for automation is available or many cells are desired, then HIVE or 10X yield good results. In general, more recently developed methods perform better. This also leads to the conclusion that older methods should be phased out, and that the development of single cell RNAseq methods is still progressing considerably.

Single-Cell RNA Sequencing of Donor-Reactive T Cells Reveals Role of Apoptosis in Donor-Specific Hyporesponsiveness of Kidney Transplant Recipients.

After kidney transplantation (KT), donor-specific hyporesponsiveness (DSH) of recipient T cells develops over time. Recently, apoptosis was identified as a possible underlying mechanism. In this study, both transcriptomic profiles and complete V(D)J variable regions of TR transcripts from individual alloreactive T cells of kidney transplant recipients were determined with single-cell RNA sequencing. Alloreactive T cells were identified by CD137 expression after stimulation of peripheral blood mononuclear cells obtained from KT recipients (N = 7) prior to and 3-5 years after transplantation with cells of their donor or a third party control. The alloreactive T cells were sorted, sequenced and the transcriptome and T cell receptor profiles were analyzed using unsupervised clustering. Alloreactive T cells retain a highly polyclonal T Cell Receptor Alpha/Beta repertoire over time. Post transplantation, donor-reactive CD4+ T cells had a specific downregulation of genes involved in T cell cytokine-mediated pathways and apoptosis. The CD8+ donor-reactive T cell profile did not change significantly over time. Single-cell expression profiling shows that activated and pro-apoptotic donor-reactive CD4+ T cell clones are preferentially lost after transplantation in stable kidney transplant recipients.

Mutation Analysis of Pancreatic Juice and Plasma for the Detection of Pancreatic Cancer.

Molecular profiling may enable earlier detection of pancreatic cancer (PC) in high-risk individuals undergoing surveillance and allow for personalization of treatment. We hypothesized that the detection rate of DNA mutations is higher in pancreatic juice (PJ) than in plasma due to its closer contact with the pancreatic ductal system, from which pancreatic cancer cells originate, and higher overall cell-free DNA (cfDNA) concentrations. In this study, we included patients with pathology-proven PC or intraductal papillary mucinous neoplasm (IPMN) with high-grade dysplasia (HGD) from two prospective clinical trials (KRASPanc and PACYFIC) for whom both PJ and plasma were available. We performed next-generation sequencing on PJ, plasma, and tissue samples and described the presence (and concordance) of mutations in these biomaterials. This study included 26 patients (25 PC and 1 IPMN with HGD), of which 7 were women (27%), with a median age of 71 years (IQR 12) and a median BMI of 23 kg/m2 (IQR 4). Ten patients with PC (40%) were (borderline) resectable at baseline. Tissue was available from six patients (resection n = 5, biopsy n = 1). A median volume of 2.9 mL plasma (IQR 1.0 mL) and 0.7 mL PJ (IQR 0.1 mL, p < 0.001) was used for DNA isolation. PJ had a higher median cfDNA concentration (2.6 ng/µL (IQR 4.2)) than plasma (0.29 ng/µL (IQR 0.40)). A total of 41 unique somatic mutations were detected: 24 mutations in plasma (2 KRAS, 15 TP53, 2 SMAD4, 3 CDKN2A 1 CTNNB1, and 1 PIK3CA), 19 in PJ (3 KRAS, 15 TP53, and 1 SMAD4), and 8 in tissue (2 KRAS, 2 CDKN2A, and 4 TP53). The mutation detection rate (and the concordance with tissue) did not differ between plasma and PJ. In conclusion, while the concentration of cfDNA was indeed higher in PJ than in plasma, the mutation detection rate was not different. A few cancer-associated genetic variants were detected in both biomaterials. Further research is needed to increase the detection rate and assess the performance and suitability of plasma and PJ for PC (early) detection.

Liquid Biopsies for Colorectal Cancer and Advanced Adenoma Screening and Surveillance: What to Measure?

Colorectal cancer (CRC) colonoscopic surveillance is effective but burdensome. Circulating tumor DNA (ctDNA) analysis has emerged as a promising, minimally invasive tool for disease detection and management. Here, we assessed which ctDNA assay might be most suitable for a ctDNA-based CRC screening/surveillance blood test. In this prospective, proof-of-concept study, patients with colonoscopies for Lynch surveillance or the National Colorectal Cancer screening program were included between 7 July 2019 and 3 June 2022. Blood was drawn, and if advanced neoplasia (adenoma with villous component, high-grade dysplasia, ≥10 mm, or CRC) was detected, it was analyzed for chromosomal copy number variations, single nucleotide variants, and genome-wide methylation (MeD-seq). Outcomes were compared with corresponding patients' tissues and the MeD-seq results of healthy blood donors. Two Lynch carriers and eight screening program patients were included: five with CRC and five with advanced adenomas. cfDNA showed copy number variations and single nucleotide variants in one patient with CRC and liver metastases. Eight patients analyzed with MeD-seq showed clustering of Lynch-associated and sporadic microsatellite instable lesions separate from microsatellite stable lesions, as did healthy blood donors. In conclusion, whereas copy number changes and single nucleotide variants were only detected in one patient, cfDNA methylation profiles could discriminate all microsatellite instable advanced neoplasia, rendering this tool particularly promising for LS surveillance. Larger studies are warranted to validate these findings.

3D chromatin reprogramming primes human memory TH2 cells for rapid recall and pathogenic dysfunction.

Memory T cells provide long-lasting defense responses through their ability to rapidly reactivate, but how they efficiently "recall" an inflammatory transcriptional program remains unclear. Here, we show that human CD4+ memory T helper 2 (TH2) cells carry a chromatin landscape synergistically reprogrammed at both one-dimensional (1D) and 3D levels to accommodate recall responses, which is absent in naive T cells. In memory TH2 cells, recall genes were epigenetically primed through the maintenance of transcription-permissive chromatin at distal (super)enhancers organized in long-range 3D chromatin hubs. Precise transcriptional control of key recall genes occurred inside dedicated topologically associating domains ("memory TADs"), in which activation-associated promoter-enhancer interactions were preformed and exploited by AP-1 transcription factors to promote rapid transcriptional induction. Resting memory TH2 cells from patients with asthma showed premature activation of primed recall circuits, linking aberrant transcriptional control of recall responses to chronic inflammation. Together, our results implicate stable multiscale reprogramming of chromatin organization as a key mechanism underlying immunological memory and dysfunction in T cells.

Machine learning-based somatic variant calling in cell-free DNA of metastatic breast cancer patients using large NGS panels.

Next generation sequencing of cell-free DNA (cfDNA) is a promising method for treatment monitoring and therapy selection in metastatic breast cancer (MBC). However, distinguishing tumor-specific variants from sequencing artefacts and germline variation with low false discovery rate is challenging when using large targeted sequencing panels covering many tumor suppressor genes. To address this, we built a machine learning model to remove false positive variant calls and augmented it with additional filters to ensure selection of tumor-derived variants. We used cfDNA of 70 MBC patients profiled with both the small targeted Oncomine breast panel (Thermofisher) and the much larger Qiaseq Human Breast Cancer Panel (Qiagen). The model was trained on the panels' common regions using Oncomine hotspot mutations as ground truth. Applied to Qiaseq data, it achieved 35% sensitivity and 36% precision, outperforming basic filtering. For 20 patients we used germline DNA to filter for somatic variants and obtained 245 variants in total, while our model found seven variants, of which six were also detected using the germline strategy. In ten tumor-free individuals, our method detected in total one (potentially germline) variant, in contrast to 521 variants detected without our model. These results indicate that our model largely detects somatic variants.

PD-L1 checkpoint blockade promotes regulatory T cell activity that underlies therapy resistance.

Despite the clinical success of immune checkpoint blockade (ICB), in certain cancer types, most patients with cancer do not respond well. Furthermore, in patients for whom ICB is initially successful, this is often short-lived because of the development of resistance to ICB. The mechanisms underlying primary or secondary ICB resistance are incompletely understood. Here, we identified preferential activation and enhanced suppressive capacity of regulatory T cells (Treg cells) in αPD-L1 therapy-resistant solid tumor-bearing mice. Treg cell depletion reversed resistance to αPD-L1 with concomitant expansion of effector T cells. Moreover, we found that tumor-infiltrating Treg cells in human patients with skin cancer, and in patients with non-small cell lung cancer, up-regulated a suppressive transcriptional gene program after ICB treatment, which correlated with lack of treatment response. αPD-1/PD-L1-induced PD-1+ Treg cell activation was also seen in peripheral blood of patients with lung cancer and mesothelioma, especially in nonresponders. Together, these data reveal that treatment with αPD-1 and αPD-L1 unleashes the immunosuppressive role of Treg cells, resulting in therapy resistance, suggesting that Treg cell targeting is an important adjunct strategy to enhance therapeutic efficacy.

Genome-wide analysis toward the epigenetic aetiology of myelodysplastic syndrome disease progression and pharmacoepigenomic basis of hypomethylating agents drug treatment response.

Myelodysplastic syndromes (MDS) consist of a group of hematological malignancies characterized by ineffective hematopoiesis, cytogenetic abnormalities, and often a high risk of transformation to acute myeloid leukemia (AML). So far, there have been only a very limited number of studies assessing the epigenetics component contributing to the pathophysiology of these disorders, but not a single study assessing this at a genome-wide level. Here, we implemented a generic high throughput epigenomics approach, using methylated DNA sequencing (MeD-seq) of LpnPI digested fragments to identify potential epigenomic targets associated with MDS subtypes. Our results highlighted that PCDHG and ZNF gene families harbor potential epigenomic targets, which have been shown to be differentially methylated in a variety of comparisons between different MDS subtypes. Specifically, CpG islands, transcription start sites and post-transcriptional start sites within ZNF124, ZNF497 and PCDHG family are differentially methylated with fold change above 3,5. Overall, these findings highlight important aspects of the epigenomic component of MDS syndromes pathogenesis and the pharmacoepigenomic basis to the hypomethylating agents drug treatment response, while this generic high throughput whole epigenome sequencing approach could be readily implemented to other genetic diseases with a strong epigenetic component.

Zeb2 DNA-Binding Sites in Neuroprogenitor Cells Reveal Autoregulation and Affirm Neurodevelopmental Defects, Including in Mowat-Wilson Syndrome.

Functional perturbation and action mechanism studies have shown that the transcription factor Zeb2 controls cell fate decisions, differentiation, and/or maturation in multiple cell lineages in embryos and after birth. In cultured embryonic stem cells (ESCs), Zeb2's mRNA/protein upregulation is necessary for the exit from primed pluripotency and for entering general and neural differentiation. We edited mouse ESCs to produce Flag-V5 epitope-tagged Zeb2 protein from one endogenous allele. Using chromatin immunoprecipitation coupled with sequencing (ChIP-seq), we mapped 2432 DNA-binding sites for this tagged Zeb2 in ESC-derived neuroprogenitor cells (NPCs). A new, major binding site maps promoter-proximal to Zeb2 itself. The homozygous deletion of this site demonstrates that autoregulation of Zeb2 is necessary to elicit the appropriate Zeb2-dependent effects in ESC-to-NPC differentiation. We have also cross-referenced all the mapped Zeb2 binding sites with previously obtained transcriptome data from Zeb2 perturbations in ESC-derived NPCs, GABAergic interneurons from the ventral forebrain of mouse embryos, and stem/progenitor cells from the post-natal ventricular-subventricular zone (V-SVZ) in mouse forebrain, respectively. Despite the different characteristics of each of these neurogenic systems, we found interesting target gene overlaps. In addition, our study also contributes to explaining developmental disorders, including Mowat-Wilson syndrome caused by ZEB2 deficiency, and also other monogenic syndromes.

Rabies virus uniquely reprograms the transcriptome of human monocyte-derived macrophages.

Macrophages are amongst the first immune cells that encounter rabies virus (RABV) at virus entry sites. Activation of macrophages is essential for the onset of a potent immune response, but insights into the effects of RABV on macrophage activation are scarce. In this study we performed high-throughput sequencing on RNA extracted from macrophages that were exposed to RABV for 48 hours, and compared their transcriptional profiles to that of non-polarized macrophages (M0), and macrophages polarized towards the canonical M1, M2a and M2c phenotypes. Our analysis revealed that RABV-stimulated macrophages show high expression of several M1, M2a and M2c signature genes. Apart from their partial resemblance to these phenotypes, unbiased clustering analysis revealed that RABV induces a unique and distinct polarization program. Closer examination revealed that RABV induced multiple pathways related to the interferon- and antiviral response, which were not induced under other classical polarization strategies. Surprisingly, our data show that RABV induces an activated rather than a fully suppressed macrophage phenotype, triggering virus-induced activation and polarization. This includes multiple genes with known antiviral (e.g. APOBEC3A, IFIT/OAS/TRIM genes), which may play a role in anti-RABV immunity.

Retrospective analysis of enhancer activity and transcriptome history.

Cell state changes in development and disease are controlled by gene regulatory networks, the dynamics of which are difficult to track in real time. In this study, we used an inducible DCM-RNA polymerase subunit b fusion protein which labels active genes and enhancers with a bacterial methylation mark that does not affect gene transcription and is propagated in S-phase. This DCM-RNA polymerase fusion protein enables transcribed genes and active enhancers to be tagged and then examined at later stages of development or differentiation. We apply this DCM-time machine (DCM-TM) technology to study intestinal homeostasis, revealing rapid and coordinated activation of enhancers and nearby genes during enterocyte differentiation. We provide new insights in absorptive-secretory lineage decision-making in intestinal stem cell (ISC) differentiation and show that ISCs retain a unique chromatin landscape required to maintain ISC identity and delineate future expression of differentiation-associated genes. DCM-TM has wide applicability in tracking cell states, providing new insights in the regulatory networks underlying cell state changes.

USP7 regulates the ncPRC1 Polycomb axis to stimulate genomic H2AK119ub1 deposition uncoupled from H3K27me3.

Ubiquitin-specific protease 7 (USP7) has been implicated in cancer progression and neurodevelopment. However, its molecular targets remain poorly characterized. We combined quantitative proteomics, transcriptomics, and epigenomics to define the core USP7 network. Our multi-omics analysis reveals USP7 as a control hub that links genome regulation, tumor suppression, and histone H2A ubiquitylation (H2AK119ub1) by noncanonical Polycomb-repressive complexes (ncPRC1s). USP7 strongly stabilizes ncPRC1.6 and, to a lesser extent, ncPRC1.1. Moreover, USP7 represses expression of AUTS2, which suppresses H2A ubiquitylation by ncPRC1.3/5. Collectively, these USP7 activities promote the genomic deposition of H2AK119ub1 by ncPRC1, especially at transcriptionally repressed loci. Notably, USP7-dependent changes in H2AK119ub1 levels are uncoupled from H3K27me3. Even complete loss of the PRC1 catalytic core and H2AK119ub1 has only a limited effect on H3K27me3. Besides defining the USP7 regulome, our results reveal that H2AK119ub1 dosage is largely disconnected from H3K27me3.

CABA-V7: a prospective biomarker selected trial of cabazitaxel treatment in AR-V7 positive prostate cancer patients.

BACKGROUND: Metastatic castration-resistant prostate cancer (mCRPC) patients with positive AR-V7 expression in their circulating tumour cells (CTCs) rarely derive benefit from abiraterone and enzalutamide. DESIGN: We performed a prospective, multicenter, single arm phase II clinical trial (CABA-V7) in mCRPC patients previously treated with docetaxel and androgen deprivation therapy. OBJECTIVE: In this trial, we investigated whether cabazitaxel treatment resulted in clinically meaningful PSA response rates in patients with positive CTC-based AR-V7 expression and collected liquid biopsies for genomic profiling. RESULTS: Cabazitaxel was found to be modestly effective, with only 12% of these patients obtaining a PSA response. Genomic profiling revealed that CTC-based AR-V7 expression was not associated with other known mCRPC-associated alterations. CTC-based AR-V7 status and dichotomised CTC counts were observed as independent prognostic markers at baseline. CONCLUSIONS: AR-V7 positivity predicted poor overall survival (OS). However, cabazitaxel-treated AR-V7 positive patients and those lacking AR-V7 positivity, who received cabazitaxel as standard of care, appeared to have similar OS. Therefore, despite the low response rate, cabazitaxel may still be an effective treatment in this poor prognosis, AR-V7 positive patient population.

Molecular and behavioral consequences of Ube3a gene overdosage in mice.

Chromosome 15q11.2-q13.1 duplication syndrome (Dup15q syndrome) is a severe neurodevelopmental disorder characterized by intellectual disability, impaired motor coordination, and autism spectrum disorder. Chromosomal multiplication of the UBE3A gene is presumed to be the primary driver of Dup15q pathophysiology, given that UBE3A exhibits maternal monoallelic expression in neurons and that maternal duplications typically yield far more severe neurodevelopmental outcomes than paternal duplications. However, studies into the pathogenic effects of UBE3A overexpression in mice have yielded conflicting results. Here, we investigated the neurodevelopmental impact of Ube3a gene overdosage using bacterial artificial chromosome-based transgenic mouse models (Ube3aOE) that recapitulate the increases in Ube3a copy number most often observed in Dup15q. In contrast to previously published Ube3a overexpression models, Ube3aOE mice were indistinguishable from wild-type controls on a number of molecular and behavioral measures, despite suffering increased mortality when challenged with seizures, a phenotype reminiscent of sudden unexpected death in epilepsy. Collectively, our data support a model wherein pathogenic synergy between UBE3A and other overexpressed 15q11.2-q13.1 genes is required for full penetrance of Dup15q syndrome phenotypes.

Effect of bovine milk fat-based infant formulae on microbiota, metabolites and stool parameters in healthy term infants in a randomized, crossover, placebo-controlled trial.

BACKGROUND: Natural enrichment of sn-2 palmitate content of infant formulae by using bovine milk fat is known to reduce formation of faecal fatty acid soaps and to improve stool consistency, but effects on gut microbiota composition are unknown. The purpose of this study was to test the influence of milk fat-based formula high in sn-2 palmitate on the infants' gut microbiota composition and to confirm the beneficial effects of the formula on formation of faecal fatty acid soaps and stool consistency. METHODS: Twenty-two healthy term, formula-fed infants were enrolled in a single-blinded randomized, crossover, placebo-controlled trial. After a 2-week run-in period, infants received either a 50% milk fat-based formula containing 39% sn-2 palmitate (MF) or a vegetable fat-based formula (VF) containing 10% sn-2 palmitate in a 2 × 2-week crossover design. Faecal microbiota composition was the primary outcome of the study. Other outcomes included faecal fatty acid soap excretion, calcium excretion, gut comfort parameters and faecal metabolites. RESULTS: Microbiota analysis showed that bifidobacteria dominated the gut microbiota of most infants. Neither alpha- nor beta-diversity was significantly influenced by the intervention. Also, abundance of metabolic pathways was independent of the intervention. The MF formula resulted in significantly lower faecal levels of palmitic acid soap (p = 0.0002) and total fatty acid soaps (p = 0.0001) than the VF formula. Additionally, calcium excretion and palmitic acid concentration were significantly (p = 0.0335) lower in stool samples after MF intervention. Furthermore, a significant physiological effect on softer stools was observed in the MF intervention compared to the VF intervention (p = 0.02). Of the 870 measured faecal metabolites, 190 were significantly different after MF and VF intervention (FDR corrected p < 0.05). Most of these were found at higher levels after MF intervention, potentially indicative of the complex structure of milk fat. Metabolites with more than twofold change between interventions were mostly lipid-derived and included several milk fat-specific fatty acids. CONCLUSIONS: Replacing part of the vegetable fat in infant formula with bovine milk fat with high sn-2 palmitate levels did not change the microbiota composition, although a reduction in faecal palmitate soaps, total fatty acid soaps and calcium excretion while improving stool consistency in the MF intervention was confirmed. In addition, 190 faecal metabolites were significantly different, many related to the fat source. TRIAL REGISTRATION: Netherlands Trial Registry Identifier: NL7815 19/06/2019.

A multiplex pedigree with pathologically confirmed multiple system atrophy and Parkinson's disease with dementia.

Multiple system atrophy is considered a sporadic disease, but neuropathologically confirmed cases with a family history of parkinsonism have been occasionally described. Here we report a North-Bavarian (colloquially, Lion's tail region) six-generation pedigree, including neuropathologically confirmed multiple system atrophy and Parkinson's disease with dementia. Between 2012 and 2020, we examined all living and consenting family members of age and calculated the risk of prodromal Parkinson's disease in those without overt parkinsonism. The index case and one paternal cousin with Parkinson's disease with dementia died at follow-up and underwent neuropathological examination. Genetic analysis was performed in both and another family member with Parkinson's disease. The index case was a female patient with cerebellar variant multiple system atrophy and a positive maternal and paternal family history for Parkinson's disease and dementia in multiple generations. The families of the index case and her spouse were genealogically related, and one of the spouse's siblings met the criteria for possible prodromal Parkinson's disease. Neuropathological examination confirmed multiple system atrophy in the index case and advanced Lewy body disease, as well as tau pathology in her cousin. A comprehensive analysis of genes known to cause hereditary forms of parkinsonism or multiple system atrophy lookalikes was unremarkable in the index case and the other two affected family members. Here, we report an extensive European pedigree with multiple system atrophy and Parkinson`s disease suggesting a complex underlying α-synucleinopathy as confirmed on neuropathological examination. The exclusion of known genetic causes of parkinsonism or multiple system atrophy lookalikes suggests that variants in additional, still unknown genes, linked to α-synucleinopathy lesions underlie such neurodegenerative clustering.